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The time-dependent thcrmo-viscoplaslie response of aerospace structures 
subjected to intense aerothermal loads is predicted using the finite-element 
method. The finite-element analysis uses the Bodner-Partom unified 
viscoplastic constitutive relations to determine rate-dependent nonlinear material 
behavior. The methodology is verified by comparison with experimental data 
and other numerical results for a uniaxially-loaded bar. The method is then used 
(1) to predict the structural response of a rectangular plate subjected to line 
heating along a centerline, and (2) to predict the thermal-structural response of a 
convcctively-cooled engine cowl leading edge subjected to aerodynamic shock* 
shock interference heating. Compared to linear elastic analysis, the viscoplastic 
analysis results in lower peak stresses and regions of plastic deformations. 


INTRODUCTION 

As hypersonic vehicles accelerate or decelerate, aerodynamic shock waves sweep across the vehicle and 
interact with local shocks and boundary layers. These shock interactions expose die vehicle surfaces to severe 
pressures and healing rates in local regions where die interactions occur. Leading edges present a particularly difficult 
design problem because of such intense, highly-localized loads. The thermal-structural response of hydrogen-cooled 
super "thermal" conducting leading edge subjected to intense aerodynamic healing was studied in [1]. A thermal- 
structural elaslo/plastic analysis with experimental verification [2) of a cowl lip design demonstrated that plastic 
effects occur and may be significant. 

Until recent years, the study of structural response at elevated temperatures due to transient thermal and 
pressure loads was not possible because of an inability to model llie viscoplastic material behavior. Traditionally, 
the thermal-structural analayis has been performed using elasto/plaslic theories. These theories lack die capabilities 
to accurately model the load rale effect, cyclic material behavior, stress relaxation and other such phenomena. Over 
the last twenty years unified viscoplastic constitutive models have evolved to meet Uicse needs. These unified 
constitutive models provide a means for analytically representing a material's response from the elastic through the 
plastic range including strain-rale dependent plastic flow, creep, s tress relaxation, and temperature effects. The first 
multi-dimensional formulation of viscoplastic constitutive model was due to Bodner and Partom [3]. Since then a 
number of constitutive models have appeared; many of these models and the related theories are summarized in 
review articles that appear in [4]. 

When unified constitutive models are implemented in finite element programs, it enables viscoplastic 
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analysis of complex structures. Finite-element methods with unified constitutive models have been under 
development for the past 15 years. A ihermo-viscoplaslic finite-element computational approach for hypersonic 
structures was presented in [5]. Analysis of convectivcly-coolcd hypersonic structures illustrated the effectiveness of 
the approach and provided insight into transient plastic structural behavior at elevated temperatures. Further 
application of the approach to high temperature structures and improvements in finite-element modeling techniques 
arc required to: (1) assist in planning thermal-structural experiments to validate the computational approach, and (2) 
uid in understanding the high-temperature bcliavior of difficult design problems such as leading edges of hypersonic 
vehicles. 

The purpose of this paper is to present the ihermo-viscoplastic analysis methodology, investigate its 
application to representative aerospace structures of current interest, and to investigate the need for such analysis. 
Tlw Bodner-Partom unified constitutive relations and the thermal -structural equations are described first. Next, the 
finite-element formulation and the computational steps are presented for the ihermo-viscoplastic analysis 
methodology. The methodlogy is evaluated by comparison of results for a uniaxially loaded bar with experimental 
data and other finite-element results. The thermo- viscoplastic response of a plate subjected to line healing along the 
centerline is then obtained. The last application example presents the transient thermal-structural response of 
convectivcly-coolcd engine cowl leading edge subjected to intense aerodynamic heating and pressure. 

BODNER-PARTOM UNIFIED VISCOPLASTIC MODEL 

The structural analysis of aerospace structures at high temperature has traditionally been carried out by 
using cither elastic or classical elaslo/plaslic analysis methodology (6J. The clastic analysis solution can t be relied 
upon once the material has yielded because clastic constitutive relations are not valid. In classical time independent 
elasto/plaslic theory, traditionally used for the post-yield analysis, the structural bcliavior is predicted using a yield 
condition (such as the von-Mises yield condition), a lime independent flow rule (such as the Prandd-Reuss equation) 
and a hardening rule. Hardening rules are used to determine reverse yielding under cyclic loading. A separate creep 
law is required to predict creep behavior. This classical elasto/plaslic theory can't predict the effect of varying strain 
rale [7] which is significant if the strain rate and temperature vary over a wide range. 

Several unified viscoplaslic models have been developed in recent years, which combine the lime 
independent and time dependent response into a single inelastic component. These models utilize internal slate 
variables as the repositories for the effect of past deformation history and do not rely on the existence of a 
conventional yield condition. A unified viscoplaslic model {3,81 is used in the present analysis because of its 
widespread use and its capability to closely simulate die lime dependent material response under various kinds of 
loading [5,7J. 

The Bodner-Partom unified constitutive model represents plastic strain rale in terms of stresses and internal 
variables. The internal slate variable rates, in turn, are given as functions of stresses, plastic work, and die 
current value of internal state variables. These equations, briefly described below, include certain physical concepts 
based on dislocation dynamics. Details of these equations arc given in [5,8] and are briefly described herein. 

The strain -displacement relation wriuen in the rale form is, 

hj = ey + ey = £ ( u,j + uj4 ) (•) 

where Eij denotes the total strain components and superscripts E and P denote clastic and plastic strain components, 
respectively. The displacement components are denoted by Ui and die dot on top of die variuble represents derivative 
with lime. The subscripts i j and j,i denote differentiation with respect to the spatial coordinates. The plastic strain 
rale is described by the evolution equation (also known as flow rule) and is given by [81, 

*5 = ilD 0 exp[-l(Z 2 /3J2) ,, l (2) 

ih 2 

where Do is the limiting shear strain rate, and n is the temperature dependent material parameter known as the strain 
rate sensitivity parameter. In the above equation, the dcviatoric stress components, Sy , are related to the stress 

components by, 
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(2a) 


ijj — Oy - ^-8ij Oki 


where Sij is the Dirac delta function. The dcviatoric stress invariant J 2 is given in terms of deviatoric stress 
components by, 


h = r Sy Sij (2b) 

The variable Z appearing in Eq. (2) is interpreted as a load history dependent parameter (known as internal slate 
variable) that represents hardness of the material with respect to plastic flow. This internal stale variable consists of 
isotropic ( Z 1 ) and directional (zP) hardening components, 

Z = Z 1 + Z D 

Tlic rate of change of the internal state variable isotropic hardening component, z\ is given by, 


Z =milZi-Z l ) Wp-AiZi 


z'-z 2 
. z, . 


r\ 


( 4 ) 


where mi is the isotropic hardening rate; Zi and Z 2 are the maximum and minimum values of Z 1 > respectively; 
constants Aiand nare temperature dependent material consutnts, known as thermal recovery parameter and thermal 
recovery exponent, respectively; and the plastic work rate,Wp , which is taken as measure of hardening, is given by, 

Wp = ofctjj (4a) 

The directional hardening component of the internal state variable, Z n , is given by, 

Z D =Pyay (5) 


where aij arc the direction cosines of the current stress stale, 

Uij = a,j / [CTklCTklll 


(5a) 


and the rale of change of the parameter Pij has die same general form as that lor isotropic hardening component (Eq. 
(4)) but has tensorial character. 


iPki Pki lU r2 Pu (5b) 

z ‘ \ Ipkjpkdi 

where m 2 is the directional hardening rale; Z 2 is the maximum value of zP* nnd A 2 and r2 are the temperature 
dependent material constants, known as thermal recovery parameter and thermal recovery exponent, respectively. 

In Eqs. (2-5), the paramclersDo, Zi, Z 3 , n, r 2 , m,, and m 2 are temperature-independent material constants; 
the parameters n, Z 2 , Zo, Aj, and A 2 are temperature-dependent material constants and all other parameters are 
problem variables. The material constants for different materials are given in [8-10] along with the detailed 
procedures for determining these constants. The particular values of these constants for the nickel-based supcralloy 
(B1900+HQ, used in the present study, are given in the Appendix. 


P,j = m 2 [Z3a l j-p,j]Wp - A2Z, 
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THERMAL-STRUCTURAL FINITE ELEMENT FORMULATION 


Governing Equations 

The thermal response of the structure is governed by the energy equation in llie form, 

pc — + qy = Q 

dt 


( 6 ) 


where p is the density; c is the specific heat; T is the temperature; Q is the internal heat generation rate; and qi.j is 
the spatial derivative of die heal flux qj in the id» direction. The conduction heal flux components, qi, are related 
to the temperature gradients by Fourier's law, 

. dT 

qi = - k ‘J — (6a) 

OXj 

where kij is the temperature-dependent thermal conductivity tensor. 

The structural response of die structure is goverend by the equilibrium equations, written in rate form [5] 
as, 


dijj + f, = 0 U) 

where f. arc die body force components per unit volume. The stress rale components arc related to the clastic strain 
rate components and die rate of change of tcnijieralure by die generalized Hooke's law, 

dy = Eijki ( eft - otki AT ) ® 

where E«u represents Hooke's tensor of elasticity parameters, and ctki are components of a tensor of thermal 
expansion coefficients. 

Finite-Element Equations and Solution Algoridim 

The flux-based finite -element approach [1] is used to derive the finite-element equations for the thermal 
analysis. The flux -based finite-element approach is used because it reduces die computational time and increases the 
solution accuracy [1,1 1]. The finite-element equations used to solve for rate of change of nodal temperatures, are in 
the form 

[Ml (T) = (R)i + (R)2 (9) 

where [M] is the capacitance matrix. The two vectors on the right hand side are associated with die fluxes within the 
element and across the element boundary, respectively. These vectors are defined by, 


(R)i 


|N1 dli (qi) 
i)xi 


(9a) 


[R)2 


[ (N) [Njds (qs) 

s 


(9b) 


where {N} are the interpolation functions; vector [q) represents die nodal heal flux in the i*h direction; (qs) vector 
represents the nodal heat flux on the boundary and includes specified heating, surface convection and surface radiation, 
symbols £1 and s on the integral sign, respectively, denote domain and boundary integrals. An explicit time 
marching scheme (with lumped capacitance matrix) is used to obtain die transient temperature response. Details of 
the matrices in Eq. (9) are given in [1]. 
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The flux-based solution algorithm is also used to derive the finite-element equations for the structural 
analysis [1 1]. The finilc-clcmcnl equations arc given by. 


(KJ{u) - {R)p+{R}r+{R}o+{R}a (1Q) 

where [K] is the element stiffness matrix. The load vectors on the right hand side of the above equation arc 
associated with the rale of change of plastic strains, prescribed temperatures, surface tractions, and body forces, 
respectively. These load vectors are defined by. 


(R)r 


{R)r 


(R)o 


Jil 


|B] t [E| {e P )dft 


[BJ t |EJ {ocJAT dft 


u 


j [NJ t (a) ds 


(10a) 

(10b) 

(10c) 


(R)b = 

where [E] is the elasticity matrix; [N] is the matrix of element interpolation functions; [B] is the strain-displacement 
matrix; ft denotes the element domain; and s denotes the boundary of the clement where tractions arc specified. 
Details of these matrices can be found in ( 12], 

To predict the structural response, the solution algorithm proceeds through the following steps: 


J U 


IN] 1 (f) dft 


00d) 


1 . At time t = 0, initialize stresses ( citj) and internal state variable components (Z 1 and Z°). 

2. Compute plastic strain rates using Eq. (2). 

3. Form finite-element matrices of Eq. (10) and assemble them to obtain global matrix equations. 

4. Solve global matrix equations [K] (u) = (F) for displacement rates. 

5. Compute strain rates from the displacement rales using Eq. (1). 

6. Compute stress rates from the generalized Hooke's law using Eq. (8). 

7. Compute the internal slate variable rates using Eqs. (4-5). 

8. Integrate stress rates and internal state variable rales to obtain the solutions of the stress components and 
internal slate variables at lime l+At. 

9. If l+Al. < tfuwi go to step 2, otherwise stop. 


The above solution algorithm has been developed into a finite-element code, which includes eight-node 
three-dimensional hexahedrai elements, four-node two-dimensional quadrilateral elements and three-node two- 
dimensional triangular elements. Ttic internal stale variable rales, plastic strain rates and other variable rate terms are 
integrated using Euler forward time integration scheme with a small lime-step. Small lime-step is required because 
the viscoplastic differential equations are very stiff in nature. The lime-step for viscoplastic analysis is chosen so 
that the computed solutions arc smooth in time and don't diverge. The time-step for the viscoplaslic analysis 
depends on the rate of applied thermo-mechanical loads. The clastic analysis solutions are computed at selected 
times, without any consideration of lime-step. 

VERIFICATION OF METHODOLOGY AND APPLICATIONS 


The application of this viscoplaslic analysis approach is presented for three different problems. The first 
application highlights simple material behavior by obtaining the thermo- viscoplaslic response of a bar at uniform 
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temperature subjected to uniaxial strain. The ihermo-viscoplaslic response of a rectangular plate under intense line 
healing is presented as the second application. A conveclivcly-cooled leading edge structure is then analyzed for 
transient thermal and viscoplaslic response. A high temperature nickel-based supcralloy (B1900+HI) is used as the 
structural material in all tire example problems. The material properties are listed in the Appendix. 

Uniaxial Loading of a Bar 

The viscoplaslic behavior of a nickel supcralloy bar subjected to different strain rates and temperatures were 
reported for the NASA Lewis HOST (HOl Section Technology) program [9,10]. The authors compared the bar's 
experimental viscoplaslic response with a finite element solution to verify the capability of die Bodner-Parlom 
constitutive model in predicting the effects of strain rale and temperatures on the stress-strain response of the 
material. In the present example, the finite-element algorithm, described in previous sections, is used to predict 
response of the bar and results are compared with numerical and experimental results from (9,10). The numerical 
results in [9,10] were obtained using the MARC finite-element code, which was modified to include a subroutine for 
Bodner-Partom unified constitutive model. The finite-clement model of the bar consists of one three-dimensional 
eight-node hexahcdral element. The bar is fixed at one end, while die oilier end is subjected to loading. 

In the first experiment, the end of the bar is subjected to tensile strain at a rale of 100 p in7in.-s at three 
different temperatures of 2060 °R. 2260 °R and 2460 °R. Computed stress-strain curves compare well with the strain 
controlled experimental data and finite element results from [9] as shown in Fig. 1. The stress-strain curve is linear 
for small values of strain, then yielding starts, and a saturation stress level is reached, above which strain increases 
without any further increase in stress. The saturation stress refers to the stress value, which corresponds to the 
material plastic flow. As the temperature of the bar is increased (see Fig. 1). the material softens and a lower value 
of saturadon stress is obtained. The elastic modulus is reduced as temperature increases, as shown in Fig. 1. 


Stress 



Fig. 1 Effect of temperarture on the stress-strain curves of a nickel-based supcralloy (B 1900+Hf) 
bar at a uniaxial strain rale of 100 p in7in.-s. 
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The effect of different applied strain rales on the bar is presented in the second experiment In this second 
experiment, the bar is kept at a constant temperature of 2060 °R. The stress-strain curves, computed using the 
present methodology, compare well with experimental data and another numerical solution [9] for three different 
strain rales as shown in Fig. 2. The material behavior is unique and linear for small strain. At higher strain levels, 
the material response is not unique at different strain rates because of the viscoplastic effect The material becomes 
harder when it is pulled at a higher strain rate (e.g. at 4000 p in./in.-s), resulting in a higher value of saturation 
stress. 

Creep phenomena is demoslratcd by applying a constant load to an isothermal bar. The strain increases 
with time without any increase in applied load, which is the basic characteristic of creep loading. Strain history 
compares well with the available experimental results and numerical solution [10] as shown in Fig. 3. The figure 
shows that an isothermal bar at 2060 °R and at a constant applied stress of 74 ksi results in a larger creep strain rate 
than the bar at 2460 °R and at a constant applied stress of 12 ksi. The effect of cyclic load at two different 
temperatures, 2060 °R and 2260 °R, is shown in Fig. 4. The bar is subjected to a maximum of ±1.5% strain at a 
rate of ±100 p in./in.-s. The virgin isothermal bar is subjected to the cyclic strain loading and two stress-strain 
curves are obtained corresponding to two different temperatures. At higher temperature, the saturation stress level is 
lower in both tension and compression. Significant plastic strain is observed at zero stress level with relatively 
more plastic strain for the bar at higher temperature. 

The above example problems highlight some of the material responses, which arc predicted using the 
Bodncr-Partom viscoplastic model applied to a uniaxially-loadcd bar. These examples also validate the present 
methodology and analysis code. 



Fig. 2 Effect of varying uniaxial strain rate on the stress-strain curves of a nickel-based superalloy 
(B1900fHf) bar at a uniform temperature of 2060 °R. 
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Fig. 3 Time dependent strain in a nickel-based supcralloy (B 1900+Hf) bar subjected to a constant 
load at two different temperatures. 



Fig. 4 Predicted stress-strain curves due to uniaxial cyclic loading of a nickel-based superalloy 
(B 1900+-HI) bar at two different temperatures and at a strain rate of 100 p in./in.-s. 
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Plate Subjected to Line Healing 

The thermo-elastic behavior of a flat plate subjected to local healing was investigated in the early 1950s 
[13]. The investigation consisted of theoretical study based on approximate analytical solutions and an experimental 
study of rectangular aluminum plates. In the experimental study, simple "tent like" steady-state temperature 
distributions were introduced by healing a rectangular plate along a centerline with a healing wire and maintaining 
constant temperatures along parallel edges by water flow through coolant tubes. Top and bottom surfaces of the 
plate were insulated to produce uniform one-dimensional, linear temperature variations between the heated centerline 
and cooled parallel edges. In the experimental study [13], in-planc plate displacements were permitted to occur freely, 
but out-of-plane displacement were constrained by restraints that forced the plate to remain flat. 

In this paper, the thcrmo-viscoplaslic behavior of a rectangular plate is studied assuming the plate is 
perfectly flat and a stale of plane stress occurs. The rectangular plate subjected to line healing along the centerline is 
shown schematically in Fig. 5. The plate is subjected to a uniform heating rale, Q, 0.73 Btu/fl-s and the line 
heating width on the plate is 0.125 in. The top and bottom surfaces of the plate are assumed to be perfectly 
insulated, and with the assumption of constant thermal properties, the transient temperature distribution cun be 
determined from a one-dimensional analytical solution [14], For this study, the temperature T(y,t) was computed 
from. 


T (y i) . 9Hnl . 19 L V _.('J£- op | - |t < 2lltl )y | i sin SatlM /.- n (») 

2k k«’ £(2 11 ,l)’ 4 / 1 2/ 

where / * 5.0 in. is the plate half-width; k is the thermal conductivity; k is the thermal diffusivity; and t is time. 
The initial temperature of the plate is assumed to be at room temperature (530 °R). Transient temperature 
distributions, calculated using eq.(l 1), are shown in Fig. 6 at four instants of time. For times up to 200 s, steep 
temperature gradients occur only over the central one-half (-2.5 < y < 2.5) of the plate. The highly localized 
temperature change is characteristic of low thermal conductivity nickel-based materials. Because of symmetry, only 
one-fourth of the plate was modeled in the structural analysis. A reference temperature of 530 °R was assumed as the 
stress-free temperature. The finite-element model and the structural boundary conditions are shown in Fig. 7, A 
stretched mesh was used with small elements along the centerline and near the edge in the region of the steep 
temperature gradients. The plate structural response was first computed assuming elastic behavior and then 
viscoplastic behavior. 
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Temperature, °R 

Fig. 6 Plate transient temperature distribution in the direction perpendicular to the centerline 


546 Nodes 
500 Quads 


B(.3, .075) 
A(.3, .025) 


Fig. 7 Finite-element model and structural boundary conditions of one -quarter of the plate 
heated along centerline at y=0. 



Temperature histories and comparisons of the elastic and viscoplastic stress (ax) histories at two locations 
(A and B as marked in Fig. 7) near the center of the plate are shown in Fig. 8. For clastic behavior, the stresses, ax, 
increase monotonically as temperature increases with time (Fig. 8). In contrast to the elastic model, the viscoplastic 
model predicts yielding early in the response with the stresses remaining nearly constant at 100 ksi for over 100 
seconds. The stress at point A decreases after time t = 120 seconds due to the higher temperature at point A and the 
softening of the material (see Appendix). A further understanding of the plate response is provided by Figs. 9-10, 
which show the variation of ax (x,0,t) and ay (7.5",y,t) for the clastic and viscoplastic analysis respectively. As 
might be expected the clastic analysis grossly over-estimates the stress levels. The viscoplastic predictions for ax 
show that the plate yields along the center line at y = 0 in a region -6 in. < x < 6 in. very early in the response at t 
= 10 s. In contrast, the viscoplastic predictions for ay along the edge at x « 7.5 in. (sec Fig. 10) show that the 
behavior is clastic early in the response, but yielding occurs perpendicular to the plate centerline at about t « 30 s. 
Contour plots of Ox stresses on the deformed shape of the plate, given in Fig. 11, show regions of maximum stress, 
ax, for elastic and viscoplastic analysis. The viscoplaslic analysis shows that the region near the centerline has 
reached saturation stress and stresses remain almost constant at about -110 ksi in this region with increasing 
temperature. 




(a) Temperature variation with lime (b) Stress variation with time 

Fig. 8 Temperature histories and comparison of clastic versus viscoplastic a* stress histories 
at two locations in a plate heated along its centerline. 


A further measure of the thermally-induced plastic behavior was derived by calculating an effective plastic 
strain, which is based on the normal and shearing plastic strain components. The effective plastic strain component 
is defined by [15], 

eSrectiv. = 2. y (ef-efr 5 + (cf-ef) 2 + (ef-ef) 2 + 3. (yf, 2 +y !' i 2 ) ( 12 ) 

Contours of the effective plastic strain are plotted at four different of times on the deformed plate (sec Fig. 12). 
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•360 ksi 


x 


-110 ksi 


(a) Elastic analysis (b) Viscoplastic analysis 

Fig. 1 1 Elastic and viscoplastic di stress contours on deformed plate (plate heated along centerline, y = 0). 





Fig. 12 Effective plastic strain contours, indicating the plastic region and the maximum plastic 
strain at different times. 
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These plots indicate the evolution of the plastically deformed regions near the centerline. The plots show that early 
in the response for 0 s < t < 50 s the largest inelastic behavior occurs near the plate's center and along the plate's 
centerline, which are primarily due to<"T x stress. However, later in the response for 150 s < t < 180 s, the largest 
plastic strain occurs on the plate centerline at the free edges at x = ±7.5 in. which is induced by the dominant (Ty 
stress at that location. These results indicate that after the plate returns to room temperature permanent contractions 
can be expected along the plate centerline at y = 0, and permanent elongations can be expected along the edges x ■ 

±7.5 in. 

Computational results for this example problem indicate that large, local temperature gradients induce 
significant plastic behavior over a narrow band along the plate centerline. The viscoplastic analysis of the plate was 
performed using a lime step of 0.1 s. The analysis required two hours of CPU lime on CRAY-2 to obtain the plate 
viscoplastic response. Computation of clastic solutions at the selected limes (one second interval) required 800 
seconds of CPU time. 

Engine Cowl Leading Edge 

Acrodynamically heated leading edges arc a critical design problem for scramjcl engine structures of 
hypersonic vehicles because of the highly localized, intense acrothcrmal loads induced by shock-shock interactions. 
The flow interaction at the scramjct engine inlet for the aerospace plane cowl leading edge is schematically shown in 
Fig. 13. The problem has been the subject of a number of studies [1,2]. A flow-thcrmal-struclural study [1] 
describes a finite-element approach used in prediction of the acrothcrmal loads, the structural temperatures, and the 
structural response of a convcclivcly-coolcd leading edge. In this study, super-thermal conductivity materials were 
evaluated, relative to their ability to alleviate elevated leading edge temperatures, temperature gradients, and attendant 
stress levels. Another study [16] evaluated an adaptive, unstructured finite-element approach for integrated fluid- 
thcrmal-structural analysis by analyzing a convcctivcly-cooled leading edge made of copper. A recent study [17] 
describes a three-dimensional thermal-structural analysis of a swept cowl leading edge subjected to interference 
heating. In the latter study [17], the computed peak stresses were beyond the elastic range, and a need for the plastic 
analysis was indicated. 
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The purpose of this example problem is to evaluate the significance of viscoplastic behavior on the thermal 
structural response of a leading edge made from the low thermal conductivity, high ductility, B1900+Hf nickel 
superalloy material. This material was selected because the material constants for viscoplastic model are available 
for a wide range of temperatures. The leading edge geometry and boundary conditions are shown schematically in 
Fig. 14. The outer edge is subject to intense, highly localized, transient aerodynamic heating and emits radiant 
energy to space. The inner surface is convcclivcly cooled by the direct impingement of a hydrogen flow at a constant 
temperature of 30 °R. The coolant convective coefficient is 7.8 Btu/ft^-s-°R. The coolant flow occurs at an internal 
pressure of 1000 psi which is assumed to be uniform over the internal surface. The leading edge has an external 
radius of 0.125 in. and a wall thickness or 0.015 in. The assumed structural boundary conditions, shown in Fig. 14, 
prohibit horizontal displacement at the supports, but permit the leading edge to expand or contract in the vertical 
direction. 

The specified aerodynamic heating on the outer surface of the leading edge is shown in Fig. 15. The heat 
flux shown is a consequence of a shock interaction that occurs as the vehicle accelerates through Mach 16. The 
interaction results from the vehicle's nose bow shock sweeping across the leading edge and intersecting the leading 
edge bow shock (see Fig. 13). The distribution shown in Fig 15a is based on experimental and computational 
studies described in [1] and the references contained therein. The assumed time history of the aerodynamic healing is 
shown in Fig. 15b. The leading edge is subjected to the intense heat flux pulses each of 6 ms duration. The shock 
on the leading is assumed stationary during this 6 ms duration, which represents the time required for the shock to 
sweep across the leading edge. In addition to the aerodynamic healing, an external surface pressure with a peak value 
of 1000 psi is also imposed. The spatial and temporal distributions of the external surface pressure are the same as 
the heat flux (see Fig. 1 5a and 1 5b). 


Aero. Heating 

390 Nodes 
626 Triangles 




0.25 


in. 


Aero, 
pressure 
(1000 psi peak) 


Internal pressure 
(1000 psi uniform) 


1 


T~ 


0.015 in. 


Fig. 14 Finite-element model of the convectively-coolcd leading edge structure with applied loads 
and boundary conditions. 
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(a) Aerodynamic healing and (b) Peak loading hisirory 

pressure profile 

Fig. 15 Aerodynamic healing and external pressure profiles along Ihe outer surface and ihe peak load history. 


The finite-element two-dimensional thermal and structural analyses were performed using the unstructured 
mesh containing 390 nodes and 626 triangles which is shown in Fig. 14. This mesh was developed in reference 116] 
using an adaptive procedure that clustered small elements in regions of steep temperature gradients. Plane strain 
condition is assumed for the structural analysis. A time step of 0.01 ms was chosen for the viscoplastic analysis. 
The computations were performed on CRAY-2 computer and the viscoplaslic analysis required 3.0 hours of epu 
time. The analysis required 1000 epu seconds to obtain elastic solutions at desired times. . . 

Temperature histories predicted at two locations on the leading edge arc shown in Fig 16. The temperature 
histories correspond to two points dirccdy under the most intense lieat flux, point A on the outer surface and point B 
the inner surface (Fig. 16a). The "saw-tooth" response at point A corresponds to the heat pulses shown in Fig. 
15b. The initial heat pulse causes the surface temperature to rise rapidly to a peak of 2550 °R at 6 ms. After the 
peak, the temperature decays to about 900 °R before the second heal pulse increases the temperature at 24 ms. The 
temperature rises again to a second peak of about 2550 °R at 30 ms. The temperature response at point B on the 
inner surface shows a relatively small increase of about 200 °R corresponding to the external heating. This small 
temperature change shows the predominant role of the internal convective cooling. Temperature contours at 6 ms are 
shown in Fig. 16b. These contours show the highly localized heat-affected region, and the very steep temperature 
drop from the outer to the inner surface. 

Using these prescribed temperatures, the structural response was computed quasi-stalically assuming: (1) 
elastic material behavior, and (2) inelastic material behavior according to die Bodncr-Partom model. 

Circumfcrcndul stress histories ut points A and B, us predicted by the clastic and viscoplastic analyses, arc 
shown in Fig. 17. At time l = 0 when die external pressure, internal pressure, and external heating arc 
instantaneously applied, local bending introduces compressive initial stress on the outer surface and tensile stress on 
die inner surface. Thereafter, dtc stress histories reflect a complex interplay between thermal stress and temperature- 
dependent inelastic material behavior. During the first heat pulse, temperature on the outer surface rises very rapidly 
introducing high compression stresses at point A. The viscoplasdc analysis predicts yielding within 1 ms. Point B 
at diis lime experiences small compressive stresses. Temperature on the outer surface rises to the peak value of 
about 2550 °R in only 6 ms during the first heat pulse. On die outer surface, the material softens with a steep 
reduction in the elastic modulus. At die same time, die steep temperature gradient introduces a significant thermal 
bending moment. After the heal pulse ends at 6 ms, die bending moment persists, causing increasing tensile 
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stresses at point B until a peak is reached at about 8 ms. The viscoplaslic analysis predictions suggest slight tensile 
yielding for point B at this time. From 6 ms until the second heat pulse occurs at 24 ms, the stress at point A 
remains roughly constant, but the stress at point B decreases as the wall temperature gradient diminishes. At 24 ms, 
the second heat pulse begins, temperature increases rapidly and the stress histories lend to repeat their first cycles. 
The axial stresses (perpendicular to analysis plane) for the elastic analysis arc much higher than the in-plane stresses 
(circumferential and radial stresses) due to plane-strain constraint, but the viscoplastic axial stresses are of comparable 
value to the viscoplastic in-plane stresses. 

During the two heat pulses, the predicted clastic and viscoplastic stresses follow each other closely. The 
most significant difference in stresses occur at 25 ms. Comparison of stress contours for the elastic and viscoplastic 
analyses at t = 25 ms is shown in Fig. 1 8. The stress contours show steep gradients from the outer to inner surface, 
and stresses are confined locally to the heat affected region as might be expected. Although the viscoplaslic and 
elastic stresses are at similar levels over much of the response, consideration of the stress levels alone does not fully 
explain the behavior. At saturation stress level, significant plastic strain occurs. This is particularly true at point A 
where the highest temperature is experienced. At very high temperatures, the saturation stress level is reduced (see 
Fig. 1) and considerable plastic strain may be caused by relatively low stresses. 

Contours of effective plastic strains at four times in the response are presented in Fig. 19 (for reference, 
yielding occurs at a strain of about 2000 p in./in. at room temperature). The figure shows that point A sustains a 
plastic strain of 12000 p inVin. over the period of 24 ms due to the first heal pulse. The second heat pulse causes 
less plastic strain due to built-in residual stresses, but a total of 17000 p in./in. accumulates after 48 ms. The 
contours suggest that the largest plastic strains occur near the outer surface although small plastic strains are 
experienced through the wall thickness. These plastic strains will result in permanent deformations and residual 
stresses after the leading edge temperatures have relumed to steady-slate conditions. 

The elastic and viscoplaslic analyses of a two-dimensional model of the leading edge provide valuable 
insight into the local thermal -structural behavior at elevated temperatures. For the nickel supcralloy material studied, 
significant permanent plastic strain results. To understand the three-dimensional viscoplaslic behavior of the leading 
edge, further study is required. 



(a) Temperature variation with lime ^ Temperature contours at 6 ms 

Fig. 16 Temperature history of two locations and the leading edge temperature contours at 6 ms. 
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Fig. 19 Effective plastic strain contours indicating the propagation of plastic region and maximum plastic strain. 


CONCLUDING REMARKS 

A thcrmo-viscoplastic computational method for structures subjected to severe local heating is used to 
investigate the transient thermal-structural behavior of several problems of current interest. The method employs the 
Bodncr-Partom unified viscoplastic constitutive model implemented in a finite-element approach for quasi-static 
thermal -structural response. 

The analysis method was verified by comparing with previous uniaxial experimental and numerical results 
using B1900fHf nickel supcralloy material. The two application problems considered were: (1) a flat plate in plane 
stress condition subjected to line heating along the plate centerline, and (2) a convcctivcly-coolcd leading edge in 
plane-strain condition subjected to localized heating representing a shock-shock interaction. Elastic and viscoplastic 
analyses were compared for each of the problems. The elastic analysis tended to over-estimate stresses and, of 
course, did not predict plastic deformations. The viscoplaslic analysis predicted local yielding and provided insight 
into the history of multi-dimensional inelastic behavior induced by transient heating. The viscoplastic analyses are 
relatively expensive in comparison to clastic analysis and lack experimental validation. Yet this study shows that 
such analyses can have a valuable role in understanding complex transient structural behavior at elevated 
temperatures. 
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APPENDIX 

Material properties for nickel-based supcralloy (B1900+HI) 


Thermal Constants 


Density (p) = 0.283 Ibm/in^ 


Thermal conductivity (k) and Specific Heat (c) 
Temperature (°R) k (Btu/in-s-°R) 


0 

500 

1500 

3000 


1.23x10-4 

1.80x10-4 

3.25x10-4 

6.00x10-4 


c (Btu/lbm^R) 
0.12 
0.135 
0.177 
0.26 


20 


1. Elastic Constants 

Temperature (°R) 

492.0 

762.0 

1032.0 

1302.0 

1572.0 

1842.0 

2112.0 

2382.0 

2652.0 

2. Thermal Expansion Coefficient (a) 

Temperature (°R) 

492.0 

852.0 

1212.0 

1572.0 

1932.0 

2292.0 

2652.0 

3. Viscoplastic constants 
Temperature-independent constants 

mi= 1.8616x10*3 psH; 
Zi = 43.51 x 10^ psi; 
n = r 2 = 2; 

Temperature-dependent constants 


Tempcraturc(°R) 

1860 

n 

1.055 

Zo=Zi (ksi) 

391.6 

A, = A 2 (s-1) 

0.0 


E (psi x J 06 ) 

G(psi x 

28.88 

12.54 

28.91 

12.22 

28.31 

11.94 

27.09 

11.62 

25.30 

11.14 

22.97 

10.40 

20.13 

9.31 

16.82 

7.76 

13.07 

5.65 


a (in/in-°R) 
0.638 x 10*5 
0.744 x 10*5 
0.755 x 10-5 
0.783 x 10*5 
0.855 x 10-5 
0.888 x 10*5 
0.922 x 10-5 


m 2 = 1.0480 x 10*2 

psi* 

i 


Z 2 = 16.67 x 104 

psi 



I> )= 10000 

s-1 



2060 


2260 

2460 

1.03 


0.85 

0.70 

348.0 


275.5 

174.0 

0.0055 


0.02 

0.25 
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